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ABSTRACT 



An interactive system for producing X-ray fluoroscopic 
images determines X-ray tube photon count and voltage 
for producing acceptable quality images while minimiz- 
ing X-ray radiation dosage to a subject. An image is 
created and a signal to noise (S/N) ratio is estimated 
from the image, assuming a Poisson model is assumed 
for the X-ray image. The S/N ratio is determined by 
solving several simultaneous equations and the photon 
count is estimated to produce an image with a desired 
S/N ratio. Subsequent X-ray fluoroscopy images are 
produced with the optimum photon count Q, thereby 
reducing X-ray dosage. The optimization is repeated 
periodically to readjust the system. 

3 QaiiDS, 14 Drawing Sheets 
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FIG. 7 
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FIG. 9 
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METHOD OF ACHIEVING REDUCED DOSE 
X-RAY FLUOROSCOPY BY EMPLOYING 
STATISnCAL ESnMAXION OF POISSON NOISE 

5 

CROSS REFERENCES TO RELATED 
APPLICATIONS 

This application is related to U.S. patent application 
X-Ray Dose Reduction in Pulsed Systems by Adaptive 
X-Ray Pulse Adjustment Ser. No. U.S, Pat. No. 
5,319,696 issued Jun. 7, 1994 "X-Ray Dose Reduction in 
Pulsed Systems by Adaptive X-Ray Pulse Adjustment" 
by A. Abdel-Malek, S. Roehm, J. Bloomer by Aiman A. 
Abdel-Malek, Steven P. Roehm and John J. Bloomer, 
X-Ray Fluoroscopy System For Reducing Dosage Em- 
ploying Iterative Power Ratio Estimation Ser. No. 
07,956,203, filed Nov. 5, 1992, now U.S. Pat No. 
5,293,415 by Richard I. Hartley. Aiman A. Abdel- 
Malek and John J. Bloomer and Reduced Dose X-ray 
Fluoroscopy System Employing Transform-based Esti- 20 
mation of Poisson Noise Ser. No. 07/971,65 Ifiled Nov. 
5, 1992, now allowed, by Richard Ian Hartley all as- 
signed to the present assignee, hereby incorporated by 
reference. 

23 

BACKGROUND OF THE INVENTION 

1. Field of the Invention. 

This invention relates to fluoroscopic imaging and 
more spedfically to reduction in patient X-ray dosage 
during imaging. ^ 

2. Description Of Related Art 

An X-ray procedure, known as fluoroscopy, creates a 
series of internal images of a subject Conventional 
pulsed systems produce each image by transmitting an 
X-ray pulse or other ionizing radiation from one side of 35 
the subject and detecting the transmitted radiation or 
shadow at an opposite side of the subject. The intensity 
of an X-ray radiation beam can be described by the 
following equation: 

40 

from p. 103 of Imaging Systems for Medical Diagnostics 
by Erich Krestel, Siemens Akdengesellschaft, Berlin 
and Munich, where E is the quantum energy of the 45 
X-ray photons, Jo(E) is the intensity at energy £ of an 
incident X-ray beam, ^(x,E) is the linear attenuation 
constant which changes along a direction of the lay x, 
and ^(x3) changes with photon energy E. 

Different tissues exhibit different linear attenuation as 50 
a function of X-ray photon energy E, thereby exhibiting 
different X-ray beam intensities J after transmission 
through the tissue. Adjusting the X-ray photon energy, 
therefore, can change the relative X-ray beam intensi- 
fies as they pass through different tissue types, leading 55 
to increased contrast in an image. 

The difference in intensity between the incident 
X-ray radiation Jo and the transmitted intensity J is 
proportional to the dose absorbed by the subject being 
imaged. Compton scattering and photoelectric absorp- 60 
tion accoimt for the majority of the energy absorbed by 
the subject in the spectrum used for conventional X-ray 
imaging as described on p. 27 of Medical Imaging Sys- 
tems by Albert Macovski, 1983 Prenttce-Hall, Engd- 
wood Cliffe. N J. 0763X 65 

In fluoroscopic systems, the radiation is pulsed at a 
rate to produce a continuous sequence of images, caus- 
ing the dosage to become quite large. Fluoroscopy is 



commonly used in order to position a catheter or similar 
invasive device inside a subject Since these procedures 
may take a long time, the acquired radiation accumu- 
lates to a large total dose. A primary goal of diagnostic 
and interventional X-ray fluoroscopic procedures is to 
provide an accurate diagnosis while reducing the dose 
received by the subject and medical staff. 

Attempts have been made to reduce dose absorbed by 
the subject and medical staff during fluoroscopic proce- 
dures. These attempts can be classified into three cate- 
gories: 

(1) mechanical redesign of elements of an X-ray sys- 
tem such as the X-ray grid, grid cover, scintiUator, 
table top, cassette front etc. to reduce scattering; 

(2) the use of protective gear (e.g., gloves and glasses, 
although the use of lead gloves hampers the ability 
to perform the fine movements necessary for cathe- 
ter placement ); and 

(3) control of X-ray tube parameters. 

The X-ray tube parameters that may be varied to 
reduce X-ray dosage include the following: 

a) the X-ray tube voltage, which affects the photon 
energy of the X-rays; 

b) the filament current I/?/, which affects the rate of 
emission of X-ray photons; 

c) the pulse duration T; and 

d) the pulse rate. 

Reduction of the filament current or the pulse dura- 
tion has the effect of decreasing the exposure in each 
frame but at the cost of diminished image quality. The 
image quality is dependent on the total photon count 
per tmit area, referred to herein as "photon count". The 
photon count is equal to the product of the photon rate 
(determined by filament current l/ji) and the pulse dura- 
tion T. 

Pulse duration T has been reduced to limit the radia- 
tion dose as described in Effect of Puked Progressive 
Fluoroscopy on Reduction of Radiation Dose in the Car- 
diac Catheterization Laboratory, by D. Holmes, M. Won- 
drow, J. Gray, R. Vetter, J. Fellows, and P. Julsrud, 
Journal American College of Cardiology, vol. 15, no. 1, 
pp. 159-162, January 1990,and hereby incorporated by 
reference. 

Imaging by reduced pulse rate has the advantage of 
maintaining the important diagnostic signal at its origi- 
nal high contrast level for a given dosage, but does not 
collect as many frames. However, the fixed rate reduc- 
tion methods produce visible jerky motion artifacts. 
These artifacts may also introduce time delays between 
a physician's actions and viewed results (e.g., moving a 
catheter or injecting radio-opaque dye). 

A technique for imaging using reduced pulse rates 
triggered by the subject's organ activity was disclosed 
in U.S. patent application "Fluoroscopic Method with 
Reduced X-Ray Dosage" Ser. No. 07/810,341 by Fathy 
F. Yassa, Aiman A. Abdel-Malek, John J. Bloomer, 
Chukka Srinivas fUed Dec. 9, 1991 assigned to the pres- 
ent assignee and hereby incorporated by reference. 
Although this technique reduces dosage by reducing 
the pulse rate, it do» not adjust the power transmitted 
by the X-ray source which may further reduce dose. 

Incorrectly reducing the power transmitted by the 
X-ray source may lead to poor quality images with 
reduced diagnostic content — the image may be charac- 
terized by global graininess and low contrast about 
important features such as the catheter, balloon, vessel 
boundaries, etc. Attempts to improve signal-to-noise 
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(S/N) ratio via noise reduction filters affect the overall FIG. 5 is a graph of signal-to-noise ratio vs. mean 

image quality by averaging-out the noise contribution pixel intensity as calculated by the first Fotirier method, 

and result in the image being of questionable value since FIG. 6 shows an original image used in simulations, 

the diagnostic information is less exact at lower doses FIG. 7 is a degraded Poisson image with signal-to- 
than at higher doses. 5 noise ratio Ps'/P/v^= — 16.8 dB created from FIG. 6. 

The X-ray tube voltage and current necessary to FIG. 8 is a degraded Poisson image with 

produce a high quality image also depend on the area of Ps'/Pyv= — 10.4 dB created from FIG. 6. 

the body under study. It is well known that different FIG. 9 is a degraded Poisson image with 

tissue types attenuate X-rays differently. For example Ps'/Piv= —7.3 dB created from FIG. 6, 
bone is quite dense, requiring high-energy X-ray pho- ^0 FIG. 10 is a degraded Poisson image with 

tons for penetration, while fat is quite transparent to Ps'/PiV= —0.3 dB created from FIG, 6. 

high-energy photons. Fat requires lower-energy X-rays FIG. 11 is a degraded Poisson image with 

to retrieve an image with good definition of the embed- Ps'/Piv=2.7 dB created from FIG. 6. 

ded features (e.g., contrast). FIG. 12 is a graph of signal-to-noise ratio vs. meam 

Since conventional fluoroscopy systems may incor- pi^cl intensity as calculated by the second transform 

rectly calculate X-ray tube voltage and photon count, method also having the graph of FIG. 5 superimposed 

subjects may be exposed to more radiation than is neces- ^P^" i^- 

sary, or the images produced may be grainy and lack FIG- 13 is a graph of signal-to-noise ratio vs. mean 

desired contrast intensity as calculated by the second statistical 

Currently, there is a need to determine the required method. 

X-ray tube voltage and photon count accurately so as to FIG- 1* is a graph obtained by supcnmposing FIG. 

produce a high quality image, while also minimizing the ^ on FIG. 13. 

X-ray dose to the subject DETAILED DESCRIPTION OF THE 

SUMMARY OF THE INVENTION 25 INVENTION 

A system for X-ray fluoroscopic imaging of a subject ^h^ X'^^V received by a subject during the 

results in quality images with reduced radiation dosage ^^^^f ^^""^ 



to a subject by first creating several X-ray images each 



with an X-ray tube voltage and a different mitial photon 3^ ^^^^^ U is the peak X-ray tube voltoge in kUovolts. lyj/ 
count Qinit, determining image mtensitics and variances ^-ray tube fUament current in mA. and T is the 
by statistical methods based upon assumptions that the duration of the X-ray pulse in seconds. The function 
images are Poisson images and formmg sunultaneous f^y, Irj, ft) one that depends on ft, the attenuation coef- 
equations relating to the mtensities and variance of the ficient, and is affected by the density and geometry of 
images. The equations are solved by simultaneous meth- 35 ^^^^^ ^eing irradiated, tube voltage, geometry of 
ods to determine constants used in predicting an optimal x-ray system, and the image detector. The peak 
image mtensity that will provide an image of a desired ^y^^ voltage determines the energy per X-ray photon. 
S/N ratio. An optimal photon count Qcpt is determined xhe number of photons emitted in a unit area is known 
based upon the image intensity, the initial photon count ^s the photon count Q. The photon count is propor- 
Qw//. and optimal image intensity, and is employed in 40 tional to the duration of the pulse T. The photon count 
creating subsequent images having a desired S/N ratio ^^y also be increased by increasing the fUament current 
and.having a reduced X-ray dose. inamanner determined by cahbration of the X-ray tube. 
OBJECTS OF THE INVENTION ^ particular, the photon count Q is an increasing func- 
tion of the X-ray tube filament current I/f/. The bright- 
It is an object of the present invention to minhniz i^ 45 ness of an image created is proportional to the total 
X-ray dose in X-ray fluoroscopic imagmg by employing photon count Q. In order to image moving structures, 
statistical methods m approximating X-ray parameters. the time of exposure may be reduced from seconds to a 
It is another object of the invention to provide Pois- few milliseconds. Therefore, the fUament current must 
son images having a desired signal-to^oise ratio using a be increased in order to produce an image of sufficient 
minimum of X-ray radiation. 50 brightness. 

BRIEF DESCRIFnON OF THE DRAWINGS „ '^^ * ''^^.'^ ^^^'^^ '° m exposure 

R and the amount of radiation absorbed by the subject 

The features of the invention beUeved to be novel are at a given X-ray energy. With exact total attenuation 
set forth with particularity in the appended claims. The and geometry unavailable, a maximum exposure Kmax* 
invention itself, however, both as to organization and 55 which corresponds to a Tnaximum photon count Qmax, 
method of operation, together with further objects, may is minimized rather than dose. The exposure is propor- 
best be understood by reference to the following de- tional to the photon count, where the constant of pro- 
scription taken in conjunction with the accompanying portionality 4>(E) depends on the photon energy, E. In 
drawing in which: other words, Q=R<KE) where 4>(E) is a function that 

FIG. 1 is a block diagram illustrating the operation of 60 may be determined through lookup tables. A graph 

a conventional X-ray system. showing a suitable function <I> is given on p. 79 of 

FIG. 2 is a graph of linear X-ray attenuation coefiici- Macovski supra. Pages 78 to 80 of Macovski provide a 

ents vs. X-ray photon energy for muscle, fat and bone. discussion of the relationship between dose, exposure R 

FIG. 9 is a block diagram of a reduced dose X-ray and photon count Q. 

system according to the present invention. 65 The X-ray tube voltage range is based on: 

FIG. 4 is a two-<limensional Fourier space represen- (1) The object to be examined; and 

tation of estimations of respective noise and signal re- (2) contrast range necessary for the diagnosis (for 

gions. example, an exposure of the "bony thorax" requires 
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66 KVp in order to diagnose the bone structure, a* field of view control unit 18, causing a field of view 
whereas 125 KVp is required if the lung structure is mask 20 to be opened, allowing X-rays from X-ray tube 
to be diagnosed). 2 to pass through subject 10 and to image plane 11. 
TTie X-ray tube voltage determines, in connection control unit 14 can be controlled to cause current 
'^^Z'^^'fy^P^^^ten.ti^^ $ source 5 to pulse the current to control voltage source 
FIG. 1. illustrates an X.«y tube comprises a cofl 3 ^ ^^^^^ ^-ray tube 2, effectively 
and a pair of plates 4<i and 4*. A current source 5 pro- • • *v v *i u w • 
vides thTfilament current, which passes through a coil P'^^^S the X-ray radiation through subject 10. The 
3 causing a number of electrons 7 to "boil-off' the coil f^^^ P^^® ^ ^ ^ averager 
3. A voltage source 6 creates a voltage difference be- 10 ^ which averages the »gnal over pulse time T for each 
tween plates 4a and 4b. Electrons 7 are repelled by point of image plane 11 and provides this signal to con- 
negatively charged plate 4a to positively charged plate t^ol '""t 14. Control unit 14 constructs an image which 
4b and accelerate at a rate proportional to the voltage is displayed on a monitor 22. 

difference applied by voltage source 6. Electrons 7 Each pixel i of the created image is composed of a 

collide with plate 4b and decelerate, causing the kinetic 15 noise component n and a signal component s, or 

energy of electrons 7 to be translated into electromag- i=s+n. The signal component s of the image represents 

netic photons 8. The energy of each photon, (proper- the ideal image that would resiilt from an elimination of 

tional to the frequency of the electromagnetic radia- all noise sources, and the noise component is equal to 

tion), is proportional to the velocity of each electron 7 j^g^ ^he difference between the actuai image and the 

as it collides with plate 4b. The frequency of the electro- 20 imaffe. 

magnetic «diation is related to its ability to penetrate ^moo s is the expected value of i. it foBows that the 

material objects. The number ofelectrons 7 that boil off «^ ♦^v ♦i^ * *i. • i t» j ■ 

coil 3 is relited to the filament current passing through ^"'f ^*~°f f ^ si^ power Vs and no«e power 

coil 3. Photons 8 emitted from plate 46 are directed ^'*%%^^fj°'^'^ ^"^^ ™ 

through a subject 10 to be imaged. Photons that pass 2S 

through subject 10 are then recorded at a recording ^ . j^Lj-^t^j* 

plane 11. Recording plane 11 may comprise photo- , "^^^ i"" '° distributed urn- 

graphic material sensitive to X-rays, or an array sensi- ^ ^^^^^^ °f frequency space, whereas the 

tive to X-rays, which is used to capture an image. ^^^nal power is concentrated m the low frequency com- 

The image captured at image plane 11 varies with the 30 ponents. 

voltage of voltage source 6 and a filament current ap- Noise is spatially uncorrelated, whereas the signal is 

phed through coil 3 from current source 5. since each correlated. The signal may be modelled as a Markov 

electron that collides with plate 4^ creates a photon that random field. The noise power can be computed from 

passes through subject 10 and illuminates a small por- the image by first performing a unitary transform on the 

tion of image plane 11. The "graininess" of the captured 35 image to create a spectrum composed of components, 

image is related to the photon count Q. The unitary transform may be a Fourier, Discrete Co- 

The difference in attenuation of photons 8 passing sine, Hadainard, Discrete Sine, Haar or Slant Trans- 
through different materials of subject 10 varies with fonn as described in "Fundamentals of Digital Image 
photon energy. This difference in attenuation between Processing" by Anil K. Jain, Prentice-Hall, (1989), As 
materials detemdnes the degree of contrast in the ere- 40 described in Jain, the above-mentioned transforms are 
ated unage. In FIG. 2 the Imeax X-ray attenuation coef- jx. ^ *„«*^- :« ir 
f ■ J. r , r ^ J i_ i*^jr umtary tr the correct scaling factor is used. If another 
ficient for muscle, fat and bone are plotted for varymg r r * - ^ *u *C * r i 
X-ray photon energy. The difference between the scahng factor is used, then ttie transforms are no longer 

curves at any given photon energy level determines the ^^^"^^ f ^ ^^^f^ P^<^.P«'- 

contrast between materials represented by the curves at 45 unchanged. The term unitary transformation 

that photon energy level. Therefore, in FIG. 1. the ^ used here to mean any unitary transform possibly 

contrast of an image acquired at image plane 11 is re- multiplied by a scaling factor. Throughout the remain- 

lated to the voltage applied across plates 4a and 4b, der of this description, a Fourier transform is described, 

The dose that subject 10 receives is related to the butanyof the above-mentioned transforms may be used 

voltage applied across plates 4a and 4b, the current SO in its place. 

passing through coil 3, and the amount of time during The high-frequency components m Fourier space 

which radiation is transmitted through subject 10. with Fybeing the ijth region, or "bin", in Fourier space 

In the X-ray system of FIG. 3 physical information of the Fourier transform of the image, are summed to 

regarding the tissue or organ of a subject 10 to be im- determine the noise power. 

aged is provided to a control unit 14 through a key- 55 piQ, 4 represents a two-dimensional Fourier space 

board 16. or pointing device 17. The operator provides ^^h the zero frequency at the center point (0,0). The 

a mimmum accepteble signal-to-noise power ratjo eross-hatched area represents a suitable choice of high 

P^AT.^ in the produced miage. Control umt 14 estab- ^ ^ ^ 

hshes mitial values for X-ray tube photon count Q/i,// ncVL-^c/ ci'>v>cii • i* *u r 

and a correspondmg fdameiit current and X-ray tube 60 ^56X256 or 512x512 pixel image the width of the 

voltage Vini, based upon conventional clinical experi- cross-hatched stnp may be chosen to equal about 32 

ence tables for this purpose. frequency bms, though this choice is somewhat arbi- 

Control unit 14 furnishes a signal to current source 5 ™s choice has been shown to give good results, 

causmg it to pass a filament current through X-ray tube however. Let the cross-hatched region be called the 

2 corresponding to the desired photon count Q. Control 65 Fourier **noise region" and be denoted Fm Let Ra7 

unit 14 also furnishes a signal to voltage source 6 cans- represent the total number of bnis in the noise region, 

ing it to produce a voltage difference across the plates Then we can compute the average noise power per bin 
of X-ray tube 2. Control unit 14 also furnishes a si^ial to by the equation: 
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underestimate of signal-to-noise ratio. This has been 
(3) observed in analysis of simulated Poisson images. 
^ ^ 2 ifz-jJ A ''^^^^ exact analysis follows. As described cor- 

\ /{Kj'>^n rectly above, the power summed over the whole of 

5 Fourier space, according to Eqs. (4) or (6) represents 
and the average total power per bin by: the contribution of signal plus noise. This may be writ- 

ten as: 

Pi^Pn-^Ps (11) 



(*) 



On the other hand, the summation over the "noise 
with the summation being over all frequency bins. region" Fivcannot be accurately equated with the noise 

By subtracting the noise power from the total power power, since there will be some signal power present in 
we obtain the signal power per bin. Finally, the signal- this region. Assume that some fraction a of the total 
to-noise ratio is equal to signal power resides in the region Fat of the frequency 

spectrum. The power estimate given by Eq. (3) is then 
Ps {Pi - Pn) <5) ^'^^ equal to noise, but to P/v+aP^. This results in 

"TJv" ^ Pn ' ^ second equation 

Often, it wiU be considered preferable to remove the ^ Pn^Pn-^^iFs (n) 
power contribution due to mean intensity, or DC . a • *u *i. r u* ^rc^ 

Jower.fromthesignalpower.stocetheDCpowerinan wb«rc P;v « the expression on *efightof^ (M 
Lage is not useful in distinguishing feature in the im- Th" ooiresponding t» Eq. (7) then ui equal 

age. This leads to a definition 



to: 



25 



(6) 



Pi-Pti (I - a)Ps (13) 



A<.k-.r^..^ ;^ ^»r,««^ If a<<l and P;v>>otPs» then Eq. 13 correctly 

As before, the ratio is aennea as 30 . • j tt tIT ^ ^ r> 

^ gives P5/P//as required. However, if Piv< <aP^, then 

• this equation reduces to (l-a)/a, the ratio of signal 

„ J (7) power outside Rn to the power in Rjv, which is not 

^ ^ what is required. 

According to Parseval's formula, the values P/ and 35 ^"^^^^^S 
P/ may be computed directly from the image, instead of 
in the Fourier domain. In particular. 



TF^Ty ~ Ps ( (I a) )"*■ ( (1 - a) )• 



(14) 



, (8) 
?/ = 2 it- 40 

The S/N ratio is assumed to be proportional to the 
^ . ^ total intensity of the image, Ps/PnsjI. The intensity is 
the sum of squares of miage mtensity values, and P /is C ^^^^ measured from the image. The left 

times the vanance of ^^^^ ^ jl4j dependent upon the inten- 

Fj^Cvaiii )=C(r2 (9) 45 ^® measured from the image. Since Pa/P^ 

^ ^ is inversely proportional to I, then: 

Some other unitary transforms, such as the Discrete 
Cosine Transform (DCT), result in a transformed re- ^ (15) 

gion defined in the first quadrant (since the DCT has no Ps ^ ^ 

negative transforms). The noise region would consist of 50 

only that part of the noise region in FIG. 4 lying in the where ro is a proportionality constant, 
first quadrant. Eq. (14) may then be written as: 

Transform-Based Noise Estimation 

The above method works very well for images that 55 Xi) = 7- + C2 
contain significant noise. For images with small 
amounts of noise, however, the noise power P^v does where 
not equal the sum over high-frequency bins 



^^^vhjF^t, i^i^l' and f(D is the value of 

since significant amounts of signal power will be present 

even in the high-frequency bins. In fact, for low noise 65 Pj- Pn 

levels, the signal power in the region F^vwill outweigh 

the contribution of noise power, causing a significant measured at the given intensity value. The values ci and 

overestimate of noise power, and hence a significant C2 are constants, independent of I. 
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Another way of writing Eq. (16) is: 



10 



(1^ 



The signal power is proportional to the square of the 
mean intensity: 



Now given two different images with different inten- 5 
sities, it is possible to solve a pair of linear equations of 
the fonn of Eq. (18) to obtain the values of C] and C2. If 
more than two images are available, then it is possible to 
solve a redundant set of eqiiations in two unknowns by 
a linear least-squares technique to get the values of ci 10 
and c2. Sections 14.0-14.3 of the book ''Numerical Reci- 
pes in C", by Wilham R Press, Brian P. Flannery, Saul 
A. Teukolsky and WilUam T. Vetterling, published by 
Cambridge University Press, Cambridge, England 
(1988)provide a good description of linear least squares 
techniques. 

Given ci and C2 and Eq. (17), a and ro may be found 
since: 



and 



the noise power is proportional to the mean intensity, 



where ki and k2 are proportionality constants. 
In addition: 



+ 2 4yc 
*7 ^ 



03) 



(24) 



I + C2 



and a — 



09) 



20 



1 +C2 



or, alternatively, if the DC power is not included in the 
signal: 



Once the values of roand a have been computed, it is 
possible to predict the value of Ps/Pat at any other 25 
image intensity according to the formula 



Ps 



.2- 



Alternatively from Eqs. (11 and 12): 



Ps 
Ptf 



Pi-Pn 
Pn^oPj 



where C is the total number of pixels in the image, ix^is 
the pixel with coordinates (x,y) and the sums are over 
all pixels in the image. The bars over P5 and P at indicate 
average power per pixel. For convenience, Eq. (25) is 
chosen to represent the total power, but the following 
discussion will apply equally well to (24) with trivial 
alterations. 
Denoting: 



35 



^ x,y ^ 



40 



as the mean pixel value, and 



^ x,y ^ 



(26) 



(27) 



as the variance of pixel values, Eqs. (22, 23 and 24) may 
be combined to obtain: 



50 



(28) 



If two independent images of an identical scene with 
the same or nearly equal intensities are available, the 
present method may still be used. The images may sim- 
ply be added together pixel by pixel to obtain an image 
of double the intensity. On an assumption that the origi- 
nal images are Poisson images, the summed image will 
be a Poisson image. 

In solving a set of redundant equations by least- 
squares techniques, it is important to realize that not all 
equations should be given the same weight For in- 
stance when solving Eqs. (16) or (IS) by least-squares 
techniques, it is important to weight the measurements 
from noisy images and relatively noise-free images ap- 
propriately, otherwise poor results may be obtained. 
The optimal weights are equal to the inverse of the 
standard deviation (square root of variance) of each 
data measurement. Since it is difficult to estimate the 
variance of the value f(I) in Eq. (16), suitable weights 
may be chosen by empirical means. It was found that 55 given intensity, T: 
the choice of equal weights w/= 1.0 in solving equations 
of the form of Eq. (16) does not give very good results. 
Far better results are achieved by choosing weights 
w/=:I, where I is the sum of pixel intensity for the i-th 
image. This is equivalent to multiplying each equation 60 
by I, thereby putting the equations in the form of Eq. 
(18). Other weights that work weU arc w,=P. 
Alternative methods for estimating the S/N ratio for 

the image are based on slightly stronger assumptions _ 

about the image. These methods are preferable, when 65 weighted byl y/^ioxJ-in. This b equivalent to 
applicable, since they do not require the computation of solving the equation 
a Fourier transform of the image. 
Two assumptions are made about the image: JciJi+Az/'iBro^J-^/z, po) 



Since T and cr^ may be computed from the image, this 
gives a single linear equation in kj and ka. Two or more 
such equations may be solved for ki and ki. Then for a 



P& 
Ps 



(29) 



In solving a redundant set of equations of the form of 
Eq. (28) appropriate weights must be chosen. It may be 
computed that the variance of <r^ as defined by (Eq. 27) 
is proportional to plus lower order terms in I. For this 
reason, each equation of the form of Eq. (28) should be 
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4) Solve Eqs. (18) and £q. (19) to obtain the value of 
. parameter ro; 

Poisson Images g^j^^ ^ ^q) to result in the image intensity re- 

For true Poisson images, the situation is even more quired to form an image with the desired signal-to- 

favorable. In fact, in this case a simple method of esti- 5 noise ratio; and 

mating the signal-to-noise ratio employs a single image. 6) recalculate the required photon count, Qopt from 

For Poissonimages, it may be shown that the average the optimal image intensity, 

noise power Pn is equal to the mean pixel value. In The second image should have signal-to-noise ratio 
other words, the constant ki in Eq. (28) is equal to 1. close to the required value and any remaining correo- 
This means that it is possible to compute Ps/Pj\rfrom a 10 tion may be computed with some accuracy, 
single image, since there is only one unknown. The Statistical Methods 

resultant formula is: 

Similar methods apply to the statistical techniques. If 
p - the image is known to be a Poisson image (a fair assump- 

=, ^ j- ^ . 15 tion for X-ray images) then Eq. (31) allows the signal- 

to-noise ratio to be computed from a single high noise 

It is possible to toke advantage of the methods used "^^^ ^^^^^^ ^^^^ 

here to make an efficioit search for the optmial photon ^ 

count at any given tube voltage U. (Sn wid © © --i 

Tlie photon count Q is proportional to the to^ 2) compute the value of I required to give the re- 

mtensity I The consent of propo^^^^^^^^^ signal-to-noise ratio fuming t^t P^/P^is 

on the voltage U and the view (of the subject). How- proporti^ to I. 

ever, once an image is taken the exact rektionship L image is truly a Poisson image, then the results 

bemg the proportionality constant of Q to I may be ^5 ^^^^^ ^^^^ However, under weaker 

detenmned This allows the photon comitQ required to ^ ji^^^. ^r more images may be employed to 

create another image with a given mtensity I to be com- produce an estimate, 

puted easUy. Second Statistical Method 

Transform Based Methods. 1) Create two of more images at different intensity 

^ ^ . - t_ 1. . 1- 30 levels 

Referrmg to FIG. 5, it can be seen that when he j) solve equations of the form of Eq. (30) (by least- 
noise level is high (of the order of 10 dB or less) the ^^^^^ j ^ ^y^^ 
estunates of S/N ratio computed using Eqs. (3), (4) and values of ki and kz 

(5) (circles) are very close to the correct values com- 3J g^i^^ ^ ^29) to compute the value of I required 

puted using the straight hnc. Thus Eqs. (3), (4) and (5) 35 achieve the desired signal-to-noise ratio, 

estimate Ps/Pjv at a relatively high noise level. This second statistical method may also be used in 

provides a first method for computing Qopt. ^ases where the image is known to be a Poisson image 

Pint: Transform Method solving Eq. (30) under the additional assumption that 

1) Create an image; k2= I. 

2) compute Ps/Pivaccording to Eqs. (3). (4) and (5) at 4Q specific choice of intensity levels may be made as in 
a low signal-to-noise ratio level; following method. 

3) compute the image intensity I required to achieve -j^^ Statistical Method 
the desired S/N ratio; and j) Create an imag^ 

4) compute required photon count. Qcpt from the 2) obtain a first estimate of Qo;,/ using steps "l"-**2" of 
optimal image mtensity. 45 the first statistical method; 

Since extrapolation over several orders of magnitude 3) fonn a second image with photon count equal to 

of image intensity may be inexact, taking two or more the estimated value of Qopi;_ 

images may provide a more exact method. 4) Compute average intensity I and variance for both 

Second Transform Method images using Eqs. (26) and (27); 

1) Create two or more images at different intensity 50 5) fonn simultaneous equations according to Eq. (30); 
leveb; 6) solve simultaneous equations to obtain the values 

2) solve equations of the form of Eq. (18) (by least of ki and kr, and 

squares method if more than two images) for ci and 7) solve Eq. (29) to compute the value of Icpt required 

C2; to achieve the desired signal-to-noise ratio. 

3) compute parameter rousing Eq. (19); 55 The resulting images will have acceptable quality and 

4) solve Eq. (20) to obtain the image intensity re- vvill be produced while minimizmg the X-ray dosage to 
quired to form an image with the desired signal-to- the subject 

noise ratio; and In order to test the methods described here, simula- 

5) calculate a required photon count, Qopt from the tion was performed on artificially degraded Poisson 
optimal image intensity. 60 images. 

A more precise choice of intensity levels is as follows: FIG. 6 shows an original image that was employed in 

Third Transform Method producing Poisson images with varying degrees of Pois- 

1) Create a first image with a givai photon count son noise- More specifically, suppose ^ represents the 
Qintt pixel value at position (x,y) in the original image. Let 

2) obtain an estimate of Qopi using steps "i*'-*»4" of 65 be a real number greater than 0. A new image was 
the first Transform method. created with pixel values ixy where for each (x,y) the 

3) form a second image with photon count equal to value ixpis a sample from a Poisson random distribution 
the estimated value of Qopt; with mean fi^ A number of such images were created 
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14 



10 



IS 



20 



with different values of 0 so that the signal-to-noise 
ratio of the resulting images ranged from about —16.8 
dB to 33 dB. 

FIO. 7 is a Poissoo image created from that of FIG. 
6, with P5'/Pjv= — 16.8 dB. This was the noisiest image 
used in the simulation. FIGS. S-11 are similar Poisson 
images with S/N ratios of - 10.4 dB, -7,3 dB, —0.3 dB 
and 2.7 dB, respectively. 

FIGS. 5, 12, 13 and 14 are plots showing the results of 
applying the various method described above to FIGS. 
8-11 and additional images. In these plots, the signal-to- 
noise ratio Ps/Pjvis computed with the DC power not 
included in the signal power. The horizontal axis repre- 
sents the mean pixel value f, measured in counts per 
pixel plotted on a logarithmic scale. The vertical axis 
represents the signal^to-noise ratio in decibels. The low- 
est quality image has approximately 10- 1 counts per 
pixel and the highest qudity image approximately 10^ 
counts per pixel. 

FIG. 5 shows the result of the first transform method. 
This shows that the estimate, although linear for high 
noise images, falls off asymptotically for low noise to a 
value (l-a)/a as predicted. 

FIG. 12 shows the result of applying the second 25 
transform method. The values of P/ and Pa^ were com- 
puted for all the images. A least squares solution to £q. 
(18) was computed and the values of ro and a were 
computed using Eq. (19). For each intensity, the value 
of Ps'/^N was computed and plotted (squares). The 30 
points are shown fitted to the line given by 
Ps'/Piv=I/ro. As is apparent the fit is very good. The 
points of FIG. 5 are also superimposed (small crosses) 
on FIG. 12 to show how they deviate from the correct 
values. 

FIG. 13 shows the results of the second statistical 
method as applied to the set of images. A least squares 
solution to £q. (30) was found. The value of k2 was 
found to equal approximately 1.001, almost exactly 
equal to the theoretical value of 1 for Poisson images. ^ 
This is not surprising, since the images were generated 
using a Poisson random number generator. The signal- 
to-noise ratio was plotted in FIG. 13 using Eq. (29) to 
compute the S/N ratio. As seen, the points fit almost 
exactly to a straight line. 

FIG. 14 shows FIG. 12 superimposed upon FIG. 13. 
This shows the difference between the results of the 
two methods. For almost all points, the values of 
Ps'/Pw^are indistinguishable as computed by the second 
transform method as compared to the second statistical 
method. 

While several presently preferred embodiments of 
the invention have been described in detail herein, many 
modifications and variations will now become apparent 
to those skilled in the art It is, therefore, to be under- 
stood that the appended claims are intended to cover all 
such modifications and variations as fall within the true 
spirit of the invention. - 
What is claimed is: 5q 
1. A method of reduced dose X-ray imaging of a 
subject comprises the steps of: 

a) selecting a minimum signal-to-noise (Ps/P^m/n) 
ratio; 

b) creating an X-ray image with an initial photon 65 
count Qinit and X-ray tube voltage U; 

c) determining an average image intensity I fi-om the 
image according to the equation: 



35 



45 



where C is the total nimiber of pixels in the image 
and ixjf is the pixels with coordkates (x,y); 
d) determining a variance <r^ from the image accord- 
ing to the equation: 



e) determining P^/Pjv from the image according to 
the following equation: 



Ps 



—1 



f) determining the value of lo;;/ required to produce an 
image with a signal-to-noise ratio at least as great as 
^s/^Nmin assuming that P^/Piv is proportional to I 
according to: 



7 - r ^-s^^^^fa 



> 



h ) determining an optimal photon count Qopt re- 
quired to produce an image with intensity Topt', and 

i) creating subsequent images with the same X-ray 
voltage U and the optimal photon count Qopt hav- 
ing a signal-to-noise ratio at least as great as 
^s/^Nmin and with a reduced X-ray dose. 

2. A method of reduced dose X-ray imaging of a 
subject comprises the steps of: 

a) selecting a minimum signal-to-noise ratio 
(Ps/PiVin/«); 

b) setting a photon count Q to a value Qm/r; 

c) creating an x-ray image with a photon count Q and 
X-ray tube voltage U; 

d) determining an average image intensity I from the 
image according to: 



SO 



where C is the total number of pixels in the image 
and ix^ is the pixels with coordinates (x,y); 
e) determining a variance ar^ of the image according 
to: 



f) creating an equation proportionality constants for 
ki and ki of the form: 

g) repeating steps "c^'-'T' for several images of dif- 
fering intensity; 

h) solving the equations of step f by simultaneous 
methods to acquire ki and ki; 
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i) determining P^/Pyv according to the following 
equation: 

j) determining an optimum image intensity Icpt re- 
quired to form an image with a signal-to-noise ratio 
greater than Ps/^Nmm according to the following 
equation: 



10 
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k) determining an optimal photon count, Qop/ neces- 
sary to produce an image with image intensity Xop^ 
based upon the^observed relationship between 
image intensity Y and photon count Q, equal to 20 
Q/n//; and 

1) creating subsequent images with the same X-ray 
voltage and the optimal photon count Qo^r having a 
S/N ratio at least as great as the PyPjvm/n and with 
a reduced X-ray dose. 25 

3. A method of reduced dose X-ray imaging of a 
subject comprises the steps of: 

a) selecting a minimum signal-to-noise (Ps/Pwmui) 
ratio; 

b) creating a first X-ray image with an initial photon 30 
count Q/n,7 and X-ray tube voltage U; 

c) determining an initial average image intensity 
i=Im//of the first image according to the equation: 



35 



where C is a number of pixels in the image,and 
is a pixel intensity at an image location x^y, with tlxe ^ 
summation covering all pixel locations; 
d) determining a variance 0*^ from the image accord- 
ing to the equation: 



45 



e) detennining PVPn from the image according to ^ 
the following equation: 



Pn 
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f) determining an initial estimated value of lopr re- 
quired to produce an image with a signal-to-noise 
ratio at least as great as P^/Pyvmin assuming that 
Ps/Pys^is proportional to I according to: 



g) determining an optimal photon count Qopi based 
upon the optimal image intensity lop;; 

h) creating a second image with the same X-ray volt- 
age U and the optimal photon count Q^opt\ 

i) determining an average image intensity 1 of the 
second image according to: 

j) determining a variance of the second image accord- 
ing to: 



k) creating an equation for proportionality constants 
ki and k2 for each of the first and second images 
according to the equation: 

I) solving the equations of step k by simultaneous 
methods to acquire ki and K2; 

m) solving, to obtain a second estimate of the opti- 
mum image intensity I^pr required to form an image 
with a signal-to-noise ratio greater than Vs/PNmint 
the equation: 



( Ps/fN, 



'fNab, \ 



n) detennining an optimal photon count, Qopt re- 
quired to produce an image with intensity lop/; and 

o) creating subsequent images with the same X-ray 
voltage and the optimal photon count Qo^r having a 
S/N ratio at least as great as Vs/^Smin and with a 
reduced X-ray dose. 



65 
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